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ABSTRACT 


The probletn of forced fluid sloshing in a partially filled spinning 
spherical tank is solved numerically using the finite element method. 

The governing equations include Coriolis acceleration, empirical fluid 
damping and spatially homogeneous vorticity first introduced by Pfeiffer. 
An exponential instability similar to flutter is detected in the present 
simulation for fill ratios below 50%. This instability appears in the 
model as a result of the homogeneous vortex assumption since the free 
slosh equations are neutrally stable in the Liapunov sense. 
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I. SUMMARY 


The Galileo dual -spin Jupiter orbiter spacecraft carries an amount 
of fluid that represents about half of the total weight. The fluid is 
contained in four tanks symmetrically located about the spin axis of the 
spacecraft. In order to evaluate the behavior of the fluid slosh in the 
spinning configuration as it affects the attitude control, an analysis 
was initiated early in the project. An analytical slosh model that 
accounts for the steady spacecraft spin was performed by Messerschmitt- 
Bdlkow-Blohm Company (MBB), at Munich, W, Germany, as part of the effort 
of manufacturing the retro-propulsion module for the Galileo spacecraft, 
The MBB analysis had several simplifying assumptions: 

(a) In order to arrive at a practically solvable set of equations, 
the vorti city was taken to be independent of the spatial 
coordinates while it is a function of time only. This model 
was first introduced by Pfeiffer at MBB in 1974. 

(b) The steady-state spin as well as angular perturbation was felt 
instantaneously in the fluid. This assumption is inconsistent 
with the mechanism of vortex transfer from the tank wall into 
the fluid. 

(c) The Coriolis acceleration was retained inside the fluid but 
was neglected in the free surface pressure boundary conditions 
in order to arrive at an exact differential form. 

The assumption of homogeneous vorti city was found to lead to fic- 
titious instabilities for certain fluid fill ratios. The inconsistency 
of this model was further substantiated by a stability analysis per- 
formed by T. K. Caughey, who showed that a Liapunov function of the free 
vibration equations had a vanishing first time derivative leading to a 
neutrally stable system. It was concluded that the MBB analysis is only 
valid for fluid fills above 55%. 

An equivalent simplified dynamic pendulum model to represent the 
fuel slosh was proposed for quick assessment of the fluid effect on the 
attitude control system. Further investigation by JPL concluded that it 
is impossible to achieve such a simplified equivalence for rotating 
tanks where Coriolis acceleration and vorticity prevail. Therefore, a 
pendulum model can only give an incomplete fluid behavior representation. 

A complete analytical investigation of the fluid behavior is now 
believed to be a very complex undertaking that for cost and schedule 
reasons is not recommended. Instead, a testing program of the complete 
spacecraft with its appropriate inertia, fluid fill, and spin rates is 
suggested to obtain an understanding of the attitude control/fluid 
interaction. 
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II. INTRODUCTION 


The problem of liquid slosh in spinning containers has lately 
received substantial attention as it is an essential step in analyti- 
cally assessing the destabilizing energy dissipation in dual-spin space- 
craft. The solution of the governing equations is involved due to the 
occurrence of the Coriolis term which destroys the self-adjointness of 
the Euler equations. An alternate set of equations that excludes the 
pressure gradient can be obtained by applying the curl operator on the 
Euler equations, thus obtaining the Helmholtz equation of vorticity. 

In 1974, Pfeiffer (Ref, 1) introduced the concept of homogenous 
vorticity; he argued that the spatially independent vorticity assumes 
the value of its average over the fluid volume. Although the assumption 
is valid for completely filled ellipsodial cavities, it is only an 
approximation for a partially filled cavity, 

The fluid perturbation velocity relative to rotating coordinates 
fixed in the center of the tank is expressed by; 


u “ -vij) - • $ + si X r 


where (j) is a velocity potential, il is the homogenous vortex vector, and 
is the Stokes-Zhukovsky potential, The process of solution starts 
with the determination of trial functions <1) satisfying the problem: 


= 0 in V 


= n ‘I’n Sp (free surface) 

f I D V] n r 


D(|) 

= 0 on S,, (wetted surface) 

0 n w 


where_o^ are the eigenvalues of the potential slosh problem. The poten- 
tial !p satisfies the problem: 

= 0 in V 


i = fxn onSp + s„ 
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Substituting the set {4>ni ^ linearized Helmholz 

equation and free surface boundary conditions and using the orthogonal- 
ity of the {$nJ set, we arrive at a system of linear ordinary differen- 
tial equations in the dependent variables 0(t) and (the general- 
ized coordinates in the trial function expansion). 

The steady spin of the spacecraft is assumed to be felt all over 
the fluid instantaneously, although the periodic angular motions from 
control and nutation are only transmitted to a very thin boundary layer 
of fluid adjacent to the wetted surface of the tank with a radial thick- 
ness on the order of /vTT, where v is the kinematic viscosity and - the 
frequency. 

The method of finite elements was adopted in deteniiining '’nl 
as well as the coefficients in the resulting ordinary differential 
equations. Evaluation of the characteristic eigenvalues indicated that 
the equations possess an instability for a range the angle r (angle 
formed by the resultant gravitational acceleration and the spin axis). 
The phenomenon is very similar to flutter in the fact that two eigen- 
values coalesce in the range of instability. As in the ease of flutter 
where bending and torsion frequencies approach each other, the present 
motion involves the free •i’Tface and vortex motions. The range of 
instability was found to iii«crease as the fluid volume ratio (Vfiuid/ 
Vtank) increased, up to ?. 5B% volume ratio, beyond which the motion was 
stable for all (i. 

A stability analysis was then performed on the original governing 
equations. A Liapunov functional was established and shown to have a 
vanishing first time derivative. This proves that the motion is neu- 
trally stable, It would, therefore, appear that the assumption of 
homogenous vorticity may be true for filling volumes more than 5B%, 
while for lesser volumes this assumption cannot possibly be true, even 
approximately. 

Great care should then be taken in using Pfeiffer's model in flow 
problems involving iritrinsic resonances. 


III. EQUATIONS OF MOTION 


Let XYZ be a rotating orthogonal coordinate system with origin at 
the CG of the spacecraft such that the Z corresponds to the spin axis. 
Let xyz be a local rotating orthogonal coordinate system with its origin 
at the center of the spinning spherical tank and parallel to XYZ. The 
angular velocity vector of XYZ with respect to inertial space is: 


“R = "o’- “ 
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where u * UqIc is the steady-state spin about Z, and "= is the perturba- 
tion angular velocity vector applied to the spacecraft from control * 
rigid body, and ela'^tic appendage motions. 

Let R be the position vector from the spacecraft cente** of mass C 
to any point P in the volume V or boundary S of the sphere. Then, 


R • Rq + r 


Rq is the radius vector from C to the sphere center o, and r is the 
radius vector from o to any point in V or on S (Figure 1). 

The total velocity vector at P(x,y,z) is: 


^ ^ (^0 ’ V “ + ”o * 


where Vq « dRo/dt, and u is the liquid particle velocity relative to 
(xyz). 

At this point it is necessary to understand the mechanisms of 
transfer of vorticity from the boundary of the tank to the contained 
fluid. Assume that the tank is motionless for t < to and that the spin 
about Z reaches its steady-state value of Qq by accelerating smoothly 
in order not to incur any ripples on the free surface. Vorticity can 
only be transmitted to the liquid by stresses tangential to the boundary 
as those generated by viscosity. It can be shown that vortex motion in 
a viscous fluid is governed by the diffusion equation. It is evident 
from the solution that the time it takes to achieve a uniform fluid 
vortex is of the order: 



where a is the radius of the sphere, and v is the kinematic viscosity of 
the liquid, furthermore, enforcing a steady periodic excitation to the 
same vorticity equation shows that only a thin boundary layer of fluid 
is entrained, having a radial thickness of: 



(4b) 

i 
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The results (4a i b) demonstrate that although a steady-state spin 
equilibrium may be reached throughout the viscous liquid after some 
transient time, small periodic rotations will never be transmitted In 
the core of the liquid and can, therefore, be neglected. 

The above argwtient explains the use of (iIq x r) rather than ("r x r) 
In equation (3), since the latter form Is inconsistent with the proper 
mechanism of vortex transmission discussed above. 

The total acceleration vector at P(x,y,z) iss 






(S) 


The centrifugal acceleration term Hq ^ ^ Include the 

perturbation angular velocity w, since this latter motion Is not felt by 
the fluid. Using the definition of the substantial derivative d/dt and 
the Identity; 


u * vu ■ ^u^ - u x (7 X u) 


in the momentum equation, wq get; 


|| + ^7u^ - u X (7 X u) + 20^ X u 

+ w X Rq + % X (Aq X R) - -^P + f 


( 6 ) 


where p and p are the fluid density and pressure, respectively, and f is 
the constant body force per unit mass. Superscript (' ■ denotes the par- , , i 

tial derivative with respect to time. For an Incompressible liquid, the ^ 

continuity equation has the form: «■ 


V . V 1 = 0 => 7 . u = 0 (7) 


We distinguish two sets of kinematic boundary conditions: 

A Neuman condition on the wetted surface enforcing 
no flow through the boundary: 


V . n = 0 ==> u • n = 0 on (8a) 
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Let F(R,t) » 0 be the equation of the free surface. Then, 


+ u • VF(R,t) *= 0 (8b) 


In addition to the kinematic relations, the dynamic condition of 
constant pressure at the free surface Sp is imposed: 


p(R,t) = constant on Sp (8c) 


In general, the relative velocity vector u can be expressed as a 
combination of an irrotational potential and a rotational term to allow 
for vorticity: 


u = 7^ + X r (9) 


where 

1^7 X V = Sp = % + ^ 

is the total vorticity. Using the definition in (10) and relation (3): 
2(S + Sq) = V X V = V X u + V X (Sp X r) 

=> 7 X u = 2(S + S^) “ 7 x(Sjj X r) = 25 

Therefore, 

-u X (v X u) = 25 X u 

Substituting (9) in (7) leads to: 


72y» =0 (10) 
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Substituting (9) in (8a), we arrive at the boundary condition for <p', 


vv» . n (n X r) . n = o 


( 11 ) 


1^- = -(r X n) • on 

The form of the boundary condition (11) suggests representing In the 
form 


<p = -(ji - n • 5i (12) 


We now introduce the concept of homogenous vortex motion (Ref. 1). 
Assume that the vorticity si is independent of the spatial coordinates 
and is only a function of time. This assumption is an approximation for 
the partially filled tank, while it is an exact description for rotating 
motions of ideal fluids totally filling spheroidal containers (Ref. 2). 
Based on this assumption, (11) can be rewritten as: 


an 


-(r X n) • 


.|i . n . |i 
an an 


on 


Since the operator on ip and associated boundary condition can be 
prescribed arbitrarily, let: 


v^ip = 

0 

in V 

. H = 
’ an 

r X n 

on S^ 


0 

in V 

. li = 

’ an 

0 

on S^ 


(13) 


The Stokes-Zhukovsky potential vector ip is used to satisfy the 
inhomogenous boundary condition (ll). The scalar potential ij' satisfies 
a homogenous Neuman condition on Sw and, as will be shown later, the 
linearized free surface condition. Substituting (13) in (10) with the 
assumption that n = n(t): 


u = -V(j> - fi • vip + n X r (14) 
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The linearized free surface condition is now derived. At first we 
define the steady-state equilibrium free surface equation by estimating 
all time-dependent quantities in the momentum equation (6): 


% X ( 5 q X R) 



+ f 


where subscript "o" denotes quantities at equilibrium. Rewriting the 
cross product in the form of a gradient and using the dynamic condition 
on the free surface (8c) : 


y - X R)^ - f • R = const. 05) 


If the body force is a thrust -go along the Z axis, then f'R = -goZk, 
where k is a unit vector along z. The steady-state free surface equa- 
tion becomes: 


K^(X,Y,Z) = + V2) - 9„Z - = 0 


Fq(X,Y,Z) = 0 describes a paraboloid of revolution, where is a con- 
stant that can be determined from the known fluid volume (Figure 2.) 

Let n be the perturbation wave height in the direction normal to 
the free surface and (?,p) its components in the Z and R directions. 
Incrementing R by p and Z by ? and retaining terms to the same order of 
smallness in (15), 


£' = (16) 

where p' is the perturbation free surface pressure. Let np and nz be 
the direction cosines of the unit normal to the free surface. Then 
(Figure 2), 


"r 


^R 


9r 9r 

= _J_^ = 9o 

9r aZ gp 


(17) 
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V 


Substituting (17) in (16) and noting that 


n = -|n„9 + 


we obtain the perturbed free surface boundary condition 




08 ) 


By definition, 


n = n • u 


on Sr 


(19) 


Substituting (14) in (19) using (13), 


n = -n • V(fi 


= M 
an 


(20) 


(20) states that the perturbation velocity normal to the free surface 
results only from the potential <j), while iji is purely geometric. Substi- 
tuting (14) in (6) using (3) and (5), 




+ ||- X r + 2fip X u = 0 


( 21 ) 


•The physical meaning of each of the terms in (21) is elaborated in 
what follows: 


at 


unsteady translational potential satisfying the 
free surface boundary condition 


J - f . R 
p 


in 2 


pressure and body forces 


convective acceleration 
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/- s X - perturbation tangential acceleration caused by 

X Kq/ r rotation about the spin axis assuming that only 

translational oscillations are transmitted to the 
liquid 



steady^state centrifugal acceleration 

tangential acceleration from vortex oscillations 

total Coriolis acceleration from spin and vortex 
motions 


In order to express the perturbation pressure_at the_free surface 
p' (equation 18) explicitly in terms of the <}), il, w, and the terms 


ft X r + 2jip X u (22) 


are omitted in (21) so that the momentum equation may be expressed as an 
exact differential. That is, 


J£’ 

p 



(w X Rq) 


r 


(23) 


The neglect of the acceleration terms (22) is consistent with the 
linearized free surface condition (18). Eliminating p'/p from (18) and 
(23), we obtain the perturbed free surface boundary condition: 


+ (r X Rq) • w (24a) 

Differentiating (24a) with respect to time and eliminating n using (20), 

h ^ ^ ‘ ^ ’ l^^o ^ ^ RqI * “ 

- (r X Rjj) . w (24b) 


An approximate solution to the slosh problem can be found by adopt- 
ing the Galerkin technique. A set of simple trial functions is sought 


TO 




that satisy the boundary conditions of the irrotational problem (equa- 
tion (24b) with si n w r: 0). A candidate set is composed of the eigen- 
functions of the problem: 



0 in V 



0 on S 

w 



1 

3t2 


0 


on Sc 
F 


(25) 


The eigenfunctions of problem (25) satisfy the orthogonality condition: 


/: 


gr, "^m *^mn 


where is the Kronecker delta. 

The functions and n in (24a) are expressed in terms of a trunca 
ted series in 


^(R) 

n(R) 



( 26 ) 


Combining (20) and (25) using (26), we obtain: 


at = y , = y !a 

an V n an ^ 9t 


\ 4> 

9 r ^n ^n 


where {Of,} are the eigenfrequencies of problem (25). 


( 27 ) 


n 


Equation (27) yields to relations between and 


^n * ’®n^^n 


Substituting (26) and (28) in (24a), 


I « 

S ^n ^nl *^n “ « ’ + (^ x R^) • w (29) 


n - 1 


Multiplying both sides of (29) by <l>f,/gD, integrating over the 
free surface, and using the orthogonality of the set. 


«n " “ S <^n*i * <=inV 

i = 1 

l^in( “ ‘^[^o]|hn( 

{^inl " ^ (®inl 
^n = /*n^9R^^ 

^n = /^i-8l^\dS = /i>,^dS 

s s 

= f (r X fi)^ dS 

r 9<l>n f ax. 

^n -■ ^Va7;^"kdS = y^^^n^dS 

= An "i 

s 

is the skew symmetric tensor form of R^. 
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The total vorticitynj ts governed by 
rotating cordi nates; ^ 


the Helmholz equation In 


dSf 

"W 



X flp ® (ijp • 
flp » ^ X V 


V) V 


( 31 ) 


(31) can be derived by applying the curl operator to both sides of the 
momentutn equation. In order to remove the spatial dependence, we define 
an average vorticity over the volume; 


flp = y dV * si + Oq 
V 


Averaging equation (31), 


+ So X Qf = I j (»p • v)v dV 

V 

The expression for 7 is found by combining (3) and (14): 


(32) 


V 


•r * % * 


Sip X 


r - S2 . Vi|; - 7(|i 


(33) 


The gradient of the first two terms in (33) vanishes identically. 
The contribution of the third terms v(3) is: 





Averaging throughout the volume and linearizing, 

■7 /<"f • 

V * 

'’ij ' F /sx7<“f • ''♦jl <''' (34) 

V ^ 

" V i <“f • “ T 

, s S ^ 


13 - 


The linearized average of the 


fourth term 



can be expressed as: 



The linearized vector equation for homogeneous vorticity is then: 


0 + (R + 





(36) 


where R is a tensor whose components are given by (34), K is the vector 
defined in (35), and is the skew symmetric tensor of vector Sq. 

The coupled equations given by (30) and (36) define the approximate 
rotational unsteady slosh motion of an ideal fluid with homogeneous 
vorticity in a spinning tank. The time history can be calculated once 
the truncated sets and are determined. 

IV. SLOSH FORCES AND MOMENTS 


The slosh force and moment vectors (F, T), acting at the spacecraft 
reference point, are given by: 


F = /ph dS = /vp dV 
S J 

f = / r X pnds = /r 


R X pn dS = J R X Vp dV 
S V 


(37) 

(38) 
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Using the expression for vp in (6) and the definitions of u in (14), 


7p 


■pjll + 1-702 + n X 0 + X u + wj^ X Rq 


+ Oq X (Oq X R) *’ f| 


u * -7<j) - (n • 7)t + 0 X r 

<}>, 


* * ■ ill 

II * 1 n 


n 


n 


■ t 


Ur ‘*’n 


Substituting (39) in (37), 

F = H(f - X (a„ X s)}+ m(r„ X j) 

- Z \\^n. 


n » 1 


5 = I J RdV « Rq + s (fluid CG) 
V 




’in “ /ir'''' ' I Vi 
V 1 s 

>1n = j *n "1 “S 


dS 


(39) 


(40) 


The terms M(f - W x (Oq x S)) represent the equilibrium constant 
thrust and centrifugarbody forces. The term M(Rq x ~) is the trans- 
lational force (constant throughout the fluid) generated by the per- 
turbation angular acceleration 6. The last bracket of (40) involves the 
slosh unsteady forces as it depends on the generalized coordinates f, 
and their derivatives: 

translation inertia of slosh 
n n 

2g X ► Coriolis force due to relative velocity 

0 n n 

d^ 5 ^a ^2 ^restoring force frcn potential energy 
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Note that the Coriolis effect is equivalent to a velocity proportional 
damping force. Substituting (39) in (38), 


T * M(1S>f + ^qSOq) + MSw 
+ (B + A)n + 2(V - u)n 

3 is the skew symmetric tensor of vector '5' in (40) 


A - 

py r r dV 
V 




B « 

pi r h f dS 
s 





p /f f dV 
V 




mu 

D - 

-P%«j3 Jf 

dV 




-/'‘(v+n*' 

V 




«n“ 

J r n dS 

S ' 

®n “ 

/: 

r 

' dS 


J h dS 

S ” 

^n 

/ 

n dS 


e„ + Ro§r 

S - RqRo 

+ 

s Ro 

u = 

[(M R„ + J) 6o]’ 

- N 6„ 5 




D - M Ro “o * 




n 

fn - “o K 

I 




(41) 


16 


Quantities superscripted by (*) are tensors, while (”) denotes the 
skew symmetric niatirlx form of a vector. The first three moment terms 
In (41) result from the first three force terms In (40), while the 
dependent contributions are born from the slosh unsteady forces. The 
^ and !i dependent terms correspond to the couples produced by the vortex 
tangential acceleration and Coriolis acceleration, respectively. 


V. PROCESS OF SOLUTION 


The eipnevalue problem (25) and the inhomogeneous potential problem 
? are solved using finite elements, Essentially, the fluid domain 
bounded by the sphere wall and the equilibrium free surface Is sub- 
divided Into Isoparametric finite elements. The Laplace equation and 
free surface boundary condition are written In a quadratic Integral 
form; 


lU) - ^ /(7*)^dV - I dS (42,) 





/{''♦, )*dV 
V 


/(r 


dS 


(42b) 


The condition that the first variation of the functional In (42a) 
be stationary leads to the Laplace equation and free surface boundary 
condition in (25). A similar process on (42b) leads to the LapUce 
equation and inhomogeneous Neumann conditions of the potential ?. 

The functional I(if>) can be discretized as; 


Kf) 



,(k)(4(k)) 


(43) 


where M* is the total number of finite elements in the fluid volume V. 
Within each element, the field variable ij> can be expressed in terms of 
its values at the nodal points using the interpolation functions 

N.(‘^^(x,y,z): 


f^‘^^(x,y,z) 


n* 


N 




(44) 
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n* denotes the number* of nodes for some element and Is the field 
variable at node 1. The condition that the functional I be stationary 

1S{ 


SI(*) - E • u 




• 0 


1 ■ l,...,n 


Using (42a) and (44), 


n* jrpN^ 0Nj 0N| aN. 3N 

j»l ax oy ay az 


f5]| 


dV 






VI, d 


which can be written In matrix form for each element as; 

£ (a<5> - b< 5))*_, >0 v(,j 




'Id 


ax ax ay ay az az/k 


over V 


q(^) s I 

\ k 


on Sr 


(45) 


where I’d denoies thej:o1 umn„vector of_the nodal values 1»'j. Furthermore, 
by splitting I* into and 1 > 2 , where ^2 contains the free surface nodes 
only, (45) becomes; 


•A, Aj- 

- X 

■Q O' 


K »2. 


0 B 

1 ^ 2 ) 


[(Bg - A|^ Ag) - XB] ^2 " 0 

i, = -a;' Aj +2 


(46) 
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Applyinci tfib above procedure to the functional with the same 

Interpolation functions N'jCxtyiz) leads to a set of inhomogeneous 
simultaneous equations of the forms 





Ai» 


(47) 


a|^^ Is identical to the matrix in (45) and 


b* » J (f X n)^ N| dS 

5 


VI. SLOSH DAMPING SIMULATION 


The energy dissipation associated with a sloshing liquid is 
governed by many complex processes that depend mainly ons 

(a) Fluid kinematic viscosity 

(b) Frequency and amplitude of excitation 

(c) Tank geometry and fluid volume 

(d) Ratio between centrifugal and thrust accelerations 

To date, no theoretical investigation exists that models even 
qualitalively the damping mechanism {Ref.3-11). It is thought in this 
work that the problem is tractable only through empiricism that relies 
largely on experiment. 

Summer and Stofan (Ref, 5) measured slosh forces on spherical tanks 
subjected to translational periodic motions of various amplitudes. The 
liquid kinematic viscosity ranged from 9.29 x 10-Z m^/sec (10"® ft^/sec) 
for water to a maximum of 9.29 x 10“4 mVsec (10”‘^ ft^/sec) for glycerine. 
For some fixed geometry and liquid properties, the nondimensional first 
mode slosh force Fg/ogD® was found to vary Tineariy with the excitation 
amplitude parameter X^/D for small Xq/D, where 

Fg ® slosh force 

p = liquid density 

g s gravitational acceleration (or thrust) 

D = tank diameter 

Xq = amplitude of periodic excitation 
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A smooth transition region follows, beyond which the slosh force 
asymptotes to a constant value depending on liquid kinematic viscosity. 

The observed trend can be explained as follows. For small amplitude 
vibrations, Fg is proportional to Xq as expected for a linear oscillator 
with constant proportional damping: 


=> F = Cp pg d3 


X 


V( 


(D 


2 _ 


1 

1)2)2 + WO) )2 


0 

D 


< 0(1) 


Cp = force coefficient function of geometry and fluid filling 
u),co^ = excitation and resonant frequencies, respectively 
? = equivalent viscous damping 


C depends only on kinematic viscosity for linear motions. As the ampli- 
tude incv'eases, free surface motions become nonlinear while the observed 
asymptotic behavior suggests that ; is proportional to X^/D (s~C^X^D); 

“ '^F 1=9“" =(t) ' 


It appears from the above limiting cases that the z variance with 
Xp/D is a hyperbola: 


O" 


5 = 


L 2 + r2/^V 

ro ^ivD/ 


(48) 


the form: 


An empirically derived expression for ? is given in (Ref. 5) in 


, . 0.131 h X 10'*\ “’9*9 , ,1. , 

'o - "zT-WTSi-l 


. > _ 1 + 0.46(2 - ha) ..... . h 

^(''a’ ' 1.46(2 ■ ' ha ' ) ^ ? 


> 1 


W = F 


for 


h < 


h is the depth of the quiescent liquid in the tank. 


( 49 ) 
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In order to determine the dependence of on the liquid and 
geometric parameters, the steady-state slosh was computed for a periodic 
excitation corresponding to the fundamental siosh resonance In transla- 
tion. A viscous damping term, proportional to the free surface normal 
velocity, was Introduced In equation (30); 


Cn + + a^2 (50) 


The damping factor ; was varied until a numerical value of the 
slosh force Fg coincided with the experimental value. A hyperbola was then 
fitted to the resulting { 5 , X^/D} set In the form: 

C . ^5/ . 10(1 . 2«„)(^r (5„ 


The frequency response spectrum was computed for a periodically 
excited spherical tank based on the damping factor expression (51). The 
geometric and liquid characterlsti 


a = 18.4 cm (7.25 In.) 

p =1000 kg/m^ (62.2 Ib/ft^) 
0.5 < n* <, 2.0 

n* is the nondimenslonal frequency 


n* = 


Figure 3 compares the present 
The correlation Is satisfactory up 
mental slosh resonance. Beyond this point, the analysis tends to under- 
estimate the slosh force, while the error increases uniformly with fre- 
quency. A possible explanation for this discrepancy is given in what 
fol 1 ows . 

Liquid slosh in spherical tanks falls in one of four different 
regimes: 

(1) At low frequencies (n* £ 1.1) the small amplitude free surface 
motion is in phase with the excitation. The liquid is mainly 
acted upon by slug forces (externally applied translational 
accelerations). The response is very similar to that of a 


s were: 



^ = 0.0083 

g = 9.81 m/sec^ (32.2 ft/sec^) 
parameter 



analysis with test data from Ref. 12. 
to n* = 1.3 just after the funda- 
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linear damped single degree of freedom oscillator having a 
mass equal to the total fluid mass. 

(2) As the frequency increases beyond n* = 1.1 > the fundamental 
resonance is approached. This results in an increase of the 
liquid mobility (degrees of freedom) and greater free surface 
amplitude. The slosh wave is broken down into small drops as 
a result of splashing against the tank wall. In addition, it 
froths due to surface cavitation and release of entrained 
vapor. These phenomena are associated with an appreciable 
increase in energy loss that is considerably higher than that 
in the linear regime. 

(3) Beyond the first resonance (n* i 1.25), the excitation and 
response are out of phase, Thi? leads to an increase in 
relative velocity between the slosh wave and the solid 
boundary. Breaking of the slosh wave as it follows the 
adverse slope of the wall occurs whenever the tank is more 
than 40% full and the wave amplitude is sufficiently large. 

It is believed that the dissipation mechanism is far more complex 
in the case of a spinning tank. Experiments were conducted at Hughes 
Aircraft in 1972 (Ref. 13) on a spinning table supported by an air 
bearing. The measured rate of energy dissipation from slosh was found 
orders of magnitude higher than that predicted analytically. In fact, 
the energy dissipation from slosh has destabilizing effects for a 
certain range of the ratio between spin and transverse inertias of the 
spacecraft. 


VII. NUMERICAL RESULTS 


Since the resultant gravitational acceleration g^ enters as a 
parameter in the definition of the nondimensional frequencies, a value 
g^ = 9.81 m/sec^ was assumed throughout the computations while the 
angle 3 : 


3 


tan 


•ll^'^ox ^ ^x^S^ 


was varied by changing both thrust g and spin rate for a mixed geom- 
etery and fluid volume. The "Galileo" tank was assumed with the follow- 
ing parameters: 

a = 0.37 m 


R„^ = 0.64 m 
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p « 872 kg/m^ (light-fluid) 
V = 0,277 X 10"® m^/sec 
'^liquid'^tank = 0.5 


The tank was spun to a constant rate (Og) on a rigid platform with 
a superimposed periodic excitation about the spin axis such that the 
equivalent maximum oscillatory amplitude at the center of the tank was 
always one percent of the tank diameter (X^/D = 0.01). Thus, 

(i) — sin tot k 


u = 0.01 ^ 

“ ''ox 


where ug and w are the excitation amplitude and frequency, respectively. 

Table 1 shows three test cases in the g range of interest. 

The variation of steady-state nondimensional circumferential force 
Fy/(p 9 gD^^o/D) with frequency parameter n* for the three cases in 
Table T are shown Figures 4,5, and 6. As expected, case A (no spin) 
exhibits one resonance at = 1.25. A similar behavior is noticed for 
case C (no thrust) except that the resonance occurs at a lower n* = 1.04. 
The response for case B is different in that three resonances are 
encountered. The first resonance is associated with a vortex dominant 
motion. The next two resonance correspond to mostly potential irrota- 
tional modes about two orthogonal axes on the free surface. A weak 
fourth resonance can be detected at n* = 1.6 and relates to the slosh 
mode with two circumferential nodes at the free surface. 

In the process of numerical experiments, a slosh instability was 
detected fluid fills ^ = (Vfiuid/Vtank) less than 0.55. The exponential 
instability was found to exist in a range of 3 angles that depends on 
the fluid fill 9. The range is largest at the lower 5 and decreases 


Table 1. Representative test cases 


Case 



g> 

m/s^ 

m/s^ 

3, 

Deg 

Kil 

A 

0.5 

2.0 

9.81 

9.81 

0.01 

B 

0.5 

31.5 

5.50 

9.81 

52 

C 

0.5 

34.0 

0.0 

9.81 

90 
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uniformly as V increases to a 9 of 0.55, beyond which the instability 
vanishes. The stability boundary (Figure 7) is independent of the 
resultant gisuvitational acceleration gj^. 

The unexpected occurrence of the instability suggested further 
investigation into its nature and mechanism. The complex characteristic 
eigenvalues of the governing system of linear differential equations 
(30) and (36) were determined for different V in the range 0 3 <. 

90 deg. The variation of the nondimensional imaginary roots Fepre'sent- 
ing the frequencies and their corresponding nondimensional real parts 
representing the danrpinn is plotted versus 3 in Figures 8 through 14 
with V as a parameter. For low V (nearly empty tank) the first fre- 
quency vanishes at 3 = 0, then rapidly increases to coalesce with the 
second resonance in the range of instability. The two frequencies then 
separate once more in the stable regime, while the lowest frequency 
approaches a finite value as 3 reaches 90 deg. The damping associated 
with the first frequency starts at zero for 3 = 0, while that corre- 
sponding to the second frequency is finite. As the instability region 
is approached, the first mode damping assumes larger negative values, 
while the second mode becomes overly damped. This behavior ceases as 
we exit the instability range. The frequency separation between the 
first two modes increases uniformly with fluid fill V, an indication of 
a weaker instability. The dashed lines represent the first two reso- 
nances of the approximate eigenvalue problem (25). The different 
negative damping parameters of the first mode are plotted against 3 
with V as a parameter in Figure 15. 

The phenomenon bears resemblance to aeroelastic flutter in which 
the bending and torsional elastic wing frequencies coalesce (Ref. 14). 

In the present model, the vortex and irrotational slosh frequencies 
couple in a manner similar to the wing frequencies in flutter. 


VIII. STABILITY ANALYSIS OF FREE SLOSHING 


A separate stability analysis on the ''free-slosh" equations was 
performed for all cavity fill ratios and all values of 0 £ 3 _< 90 deg. 

Consider Euler's equation for fluid motion in a rotating system 
(Eq. 6 with s = u = 0) with homogeneous boundary conditions given by 
(8a,b,c). The steady-state solution is given by (15). 

Let v(x,t) be the infinitesimal perturbatlonal velocity of the 
fluid, and pi(x,t) be the infinitesimal perturbational pressure. Then, 
using (6) and Tl5), the equation of perturbed motion of the fluid is: 


3V 

9t 


+ 2n X V = -— 




V . V = 0 


( 52 ) 


in V, 
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( 53 ) 


Outside of Vjj there is a "gravitational field" given by: 


ig = 0)2 [xl yj] « gic 


As a result of the perturbed motion, the free surface is perturbed. 
Let the elevation of the perturbed free surface relative to the equili- 
brium free surface Sp be denoted by c, measured normal to Sp . 

’*0 ' ^0 

From (15) and (53) it is easily seen that g is normal to the free 
surface, and is directed into the fluid. 

As before. 


V • n = 0 


on S 


w. 


The pressure on the equilibrium free surface Sp is given by: 

^0 

p, = p 9g • fi ? on Sp (correct to 0(c)) 
®Fo ° 


Consider the Liapunov functional (Ref. 15): 


1 


/ p V . V 


1 


V dv + ■S' p g 
0 2 ^ ^e 


y d 3p 

c r , 


(54) 


9e = 9e • n 


Consider the time derivative of L: 


L “ /pv. #dV„ + / 
Using (52) in (55): 


y St “’0 ' ^ P9e ?3t ^Sp 

Vo Sp^ 


(55) 


L = 


-I 


p(2s^xv’V+-^v 


/ 


vpj)dV^ + J pg,?|f ds 

Or- 0 

^0 
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Figure 3. Comparison for translational slosh 


X 

= 34 .72 NEWTONS 



Figure 4. Nondimensional tangential slosh force for case A 
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Figure 7. Stability boundary 



























